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Abstract 

The multiplicity distributions of hadrons produced in central nucleus-nucleus collisions are studied 
within the hadron-resonance gas model in the large volume limit. The microscopic correlator method 
is used to enforce conservation of three charges - baryon number, electric charge, and strangeness - in 
the canonical ensemble. In addition, in the micro-canonical ensemble energy conservation is included. 
An analytical method is used to account for resonance decays. The multiplicity distributions and the 
scaled variances for negatively, positively, and all charged hadrons are calculated along the chemical 
freeze-out line of central Pb-|-Pb (Au-|-Au) collisions from SIS to LHC energies. Predictions obtained 
within different statistical ensembles are compared with the preliminary NA49 experimental results 
on central Pb-|-Pb collisions in the SPS energy range. The measured fluctuations are significantly 
narrower than the Poisson ones and clearly favor expectations for the micro-canonical ensemble. Thus 
this is a first observation of the recently predicted suppression of the multiplicity fluctuations in 
relativistic gases in the thermodynamical limit due to conservation laws. 

PACS numbers: 24.10.Pa, 24.60.Ky, 25.75.-q 
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I. INTRODUCTION 

For more than 50 years statistical models of strong interactions have served as an 

important tool to investigate high energy nuclear collisions. The main subject of the past study 
has been the mean multiplicity of produced hadrons (see e.g. Refs. Q,fly,3). Only recently, 
due to a rapid development of experimental techniques, first measurements of fluctuations of 

formed. The growing interest in the 



10| ) is motivated by expectations of 



particle multiplicity [8| and transverse momenta [9|] were per 
study of fluctuations in strong interactions (see e.g., reviews J 
anomalies in the vicinity of the onset of deconfinement llj and in the case when the expanding 
system goes through the transition line between the quark-gluon plasma and the hadron gas [12] . 
In particular, a critical point of strongly interacting matter may be signaled by a characteristic 
power-law pattern in fluctuations 13|. Apart from being an important tool in an effort to 
study the critical behavior, the study of fluctuations in the statistical hadronization model 
constitutes a further test of its validity. In this paper we make, for the first time, predictions 
for the multiplicity fluctuations in central collisions of heavy nuclei calculated within the micro- 
canonical formulation of the hadron-resonance gas model. Fluctuations are quantified by the 
ratio of the variance of the multiplicity distribution and its mean value, the so-called scaled 
variance. The model calculations are compared with the corresponding preliminary results [14] 
of NA49 on central Pb+Pb collisions at the CERN SPS energies. 

There is a qualitative difference in the properties of the mean multiplicity and the scaled 
variance of multiplicity distribution in statistical models. In the case of the mean multiplicity 
results obtained with the grand canonical ensemble (GCE), canonical ensemble (CE), and 
micro-canonical ensemble (MCE) approach each other in the large volume limit. One refers 
lere to the thermodynamical equivalence of the statistical ensembles. It was recently found 
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161 ] that corresponding results for the scaled variance are different in different ensembles. 



and thus the scaled variance is sensitive to conservation laws obeyed by a statistical system. 
The differences are preserved in the thermodynamic limit. 

The paper is organized as follows. In Section II the microscopic correlators for a relativistic 
quantum gas are calculated in the MCE in the thermodynamical limit. This allows to take into 
account conservation of baryon number, electric charge, and strangeness in the CE formulation 
and, additionally, energy conservation in the MCE. In Section III the relevant formulas for 
the scaled variance of multiplicity fluctuations are presented for different statistical ensembles 
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within the hadron-resonance gas model. The scaled variance of negative, positive and all 
charged hadrons is then calculated along the chemical freeze-out line in the temperature- 
baryon chemical potential plane. The fluctuations of hadron multiplicities in central Pb+Pb 
(Au+Au) collisions are presented for different collision energies from SIS to LHC The results 
for the GCE, CE, and MCE are compared. In Section IV the statistical model predictions for 
the scaled variances and multiplicity distributions of negatively and positively charged hadrons 
are compared with the preliminary NA49 data of central Pb+Pb collisions in the SPS energy 
range. A summary, presented in Section V, closes the paper. New features of resonance decays 
within the MCE are discussed in Appendix A, and the acceptance procedure for all charged 
hadrons is considered in Appendix B. 



II. MULTIPLICITY FLUCTUATIONS IN STATISTICAL MODELS 

The mean multiplicities of positively, negatively and all charged particles are defined as: 

i,qi<0 i,qi>0 «,<?i7^0 

where the average final state (after resonance decays) multiplicities (Aj) are equal to: 

(AT,) = {N:) + Y,{NR){n.)R. (2) 

R 

In Eq. ([2]), A^* denotes the number of stable primary hadrons of species i, the summation 
runs over all types of resonances R, and {ni)^ = J2r ^r"'^fr is the average over resonance decay 
channels. The parameters are the branching ratios of the r-th branches, nf^ is the number of 
particles of species i produced in resonance R decays via a decay mode r. The index r runs over 
all decay channels of a resonance i?, with the requirement Ylr^r^ ~ GCE formulation 

of the hadron-resonance gas model the mean number of stable primary particles, (N*), and the 
mean number of resonances, (Nr), can be calculated as: 

where V is the system volume and gj is the degeneracy factor of particle of the species j (number 
of spin states). In the thermodynamic limit, V — » oo, the sum over the momentum states can 
be substituted by a momentum integral. The (ripj) denotes the mean occupation number of a 
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single quantum state labelled by the momentum vector p , 



^^-•^■^ exp [(ep, - /.,) IT] - a, ' 



where T is the system temperature, rrij is the mass of a particle j, e^j = y + "rn^j is a single 
particle energy. A value of aj depends on quantum statistics, it is +1 for bosons and —1 for 
fermions, while Uj = gives the Boltzmann approximation. The chemical potential Hj of a 
species j equals to: 

f^j = 1j f^Q + f^B + Sj fis , (5) 

where the particle electric charge, baryon number, and strangeness, respectively, 

while fig, hb, fJ'S are the corresponding chemical potentials which regulate the average values of 
these global conserved charges in the GCE. Eqs. ([MS]) are valid in the GCE. In the limit V ^ oo , 
Eq. dMSD are also valid for the CE and MCE, if the energy density and conserved charge 
densities are the same in all three ensembles. This is usually referred to as the thermodynamical 
equivalence of all statistical ensembles. However, the thermodynamical equivalence does not 
apply to fluctuations. 

In statistical models a natural measure of multiplicity fluctuations is the scaled variance 
of the multiplicity distribution. For negatively, positively, and all charged particles the scaled 
variances read: 

_ ((AA^_)'^> ^ ((AA^^,)') {(AN^f) 
The variances in Eq. ([6]) can be presented as a sum of the correlators: 

qi<0,qj<0 qi>0,qj>0 

{{AN,,)') = Y i^N.AN,) , (7) 

where AA^j = Aj — (Aj). The correlators in Eq. ([7]) include both the correlations between 
primordial hadrons and those of final state hadrons due to the resonance decays (resonance 
decays obey charge as well as ene.gy_tu„> conse^t.on) 
In the GCE the final state correlators can be calculated as |18l |: 

(AAT, AAr,),.,.e. = {AN*AN;)g,,,,, + [(^^r) {^^)R K')/? + i^n) {An,An,)R] , (8) 

R 
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where (Arij Anj)ji = J2r^?''^i^r''^f,r ~ {T^i)R{^j)R ■ The occupation numbers, ripj, of single 
quantum states (with fixed projection of particle spin) are equal to Upj = 0, 1, . . . , oo for bosons 
and Up J = 0, 1 for fermions. Their average values are given by Eq. (jlj), and their fiuctuations 
read: 

( (^^pj)^ ) = - ('^P.i))^) = (^P,i) (1 + "i (^P,i)) = ^p,j ■ (9) 

It is convenient to introduce a microscopic correlator, (AnpjA?2k,i), which in the GCE has a 
simple form: 

(Anpj Ank,j)g.c.e. = vlj6ij6pk. (10) 

Hence there are no correlations between different particle species, i ^ j, and/or between dif- 
ferent momentum states, p 7^ k. Only the Bose enhancement, Vpj > (ripj) for aj = 1, and the 
Fermi suppression, Vpj < (upj) for aj = —1, exist for fiuctuations of primary particles in the 
GCE. The correlator in Eq. ([8]) can be presented in terms of microscopic correlators ( ITOi) : 

{AN* AN* = Yl (^^P.J- ^^K^)9■c.e. = Yl ^Jj " 

p,k p 

In the case oi i = j the above equation gives the scaled variance of primordial particles (before 
resonance decays) in the GCE. 

In the MCE, the energy and conserved charges are fixed exactly for each microscopic state 
of the system. This leads to two modifications in a comparison with the GCE. First, the 
additional terms appear for the primordial microscopic correlators in the MCE. They refiect 
the (anti) correlations between different particles, i ^ j-, and different momentum levels, p 7^ k, 
due to charge and energy conservation in the MCE, 



(Ar2pjAnk,j)r 



6ij 6pk 



\A\ 



qiqjMqq + bibjMbb + SiSjMs. 



+ {QiSj + qjSi) Mqs - {qibj + qjbi) Mgb - {ksj + bjSi) M^s 

+ CpiCkiMe, - {qitpj + g^eki) M^, + {htpj + ft^eki) Mb, - (sjCpj + s^eki) M,, ] , 

where |A| is the determinant and Mjj are the minors of the following matrix, 

/A(g2) A(6g) A{sq) A(eg)\ 

A{qb) A{9) A{sb) A{eb) 

A{qs) A{bs) A(s2) A(es) 

\A{qe) A{be) A{se) A{e^) J 



(12) 



A 



(13) 
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with the elements, A{q^) = q^v^j , A{qb) = J^pj Qj^j^h > H<i^) = Epj Qj^pj^lj , etc. 
The sum, , means integration over momentum p, and the summation over all hadron- 

resonance species j contained in the model. The first term in the r.h.s. of Eq. f|T2l) corresponds 
to the microscopic correlator (fTOl) in the GCE. Note that a presence of the terms containing 



a single particle energy, epj = y + Tn'j, in Eq. f|T2|) is a consequence of energy conservation. 
In the CE, only charges are conserved, thus the terms containing epj in Eq. (1121) are absent. 



The A in Eq. ( fT3l) becomes then the 3x3 matrix (see Ref. [16|). An important property of 
the microscopic correlator method is that the particle number fluctuations and the correlations 
in the MCE or CE, although being different from those in the GCE, are expressed by quanti- 
ties calculated within the GCE. The microscopic correlator f|T2|) can be used to calculate the 
primordial particle correlator in the MCE (or in the CE): 

{ANi AN, = J2 (^^P.^ A^k,i)m.c.e. ■ (14) 

p,k 

A second feature of the MCE (or CE) is the modification of the resonance decay contribution 
to the fluctuations in comparison to the GCE ([8]). In the MCE it reads: 

(AN.AN,)^,,,,, = {AN*AN*)m.c.e. + ^^(iVi?) {An,An,)R + ^^(AiV* AiVR)^.,^. {n,)n 

R R 

+ J2^AN* ANn)m.c.e. {ni)R + ^{ANr ANR,)m.c.e. {n,)R {nj)R' . (15) 

R R,R' 

Additional terms in Eq. ( ITSl) compared to Eq. ([8]) are due to the correlations (for primordial 
particles) induce d by energy and charge conservations in the MCE. The Eq. f|T5l) has the same 
form in the CE [l6| and MCE, the difference between these two ensembles appears because 
of different microscopic correlators f|T2l) . The microscopic correlators of the MCE together 
with Eq. ( fT4l) should be used to calculate the correlators {AN* AN*) m.c.e. , {AN* ANR)m.c.e. , 
{AN* ANR)m.c.e. , {ANr an r,) m.c.e. entering in Eq. ^ . 

The microscopic correlators and the scaled variance are connected with the width of the 



multiplicity distribution. It can be shown 17|] that in statistical models the form of the mul- 



tiplicity distribution derived within any ensemble (e.g. GCE, CE and MCE) approaches the 
Gauss distribution: 



Pg{N) = , exp 
v/27r uj {N) 



2 UJ {N) 
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in the large volume limit i.e. (A^) oo. The width of this Gaussian, a = \/ uj {N), is 
determined by the choice of the statistical ensemble, while from the thermodynamic equivalence 
of the statistical ensembles it follows that the expectation value (N) remains the same. 

III. MULTIPLICITY FLUCTUATIONS AT CHEMICAL FREEZE-OUT 

In this section we present the results of the hadron-resonance gas for the scaled variances 
in the GCE, CE and MCE along the chemical freeze-out line in central Pb+Pb (Au+Au) 
collisions for the whole energy range from SIS to LHC Mean hadron multiplicities in heavy 
ion collisions at high energies can be approximately fitted by the GCE hadron-resonance gas 
model. The fit parameters are the volume V, temperature T, chemical potential fiB, and the 
strangeness saturation parameter 75. The latter allows for non-equilibrium strange hadron 
yields. A recent discussion of system size and energy dependence of freeze-out parameters 
and comparison of freeze-out criteria can be found in Refs. . There are several programs 
designed for the analysis of particle multiplicities in relativistic heavy-ion collisions within the 
hadron-resonance gas model, see e.g., SHARE [IS], THERMUS fl, and THERMINATOR |21|. 
The set of model parameters, V,T, fis, and 75, corresponds to the chemical freeze-out conditions 
in heavy-ion collisions. The numerical values and evolution of the model parameters with the 
collision energy are taken from previous analysis of multiplicities data. The dependence of /i^ 
on the collision energy is parameterized as [3]: /i^ {^/snn) = 1.308 GeV- (1 + 0.273 ^/SNN)^^ , 
where the cm. nucleon-nucleon collision energy, y/s^N, is taken in GeV units. The system 
is assumed to be net strangeness free, S = 0, and to have the charge to baryon ratio of the 
initial colliding nuclei, Q/B = 0.4. These two conditions define the system strange, /is', and 
electric, /ig, chemical potentials. For the chemical freeze-out condition we chose the average 
energy per particle, {E)/{N) = 1 GeV [2^- Finally, the strangeness saturation factor, 75, is 
parametrized 6j, 75 = 1 — 0.396 exp (— 1.23 T/fis) ■ This determines all parameters of the 
model. In this paper an extended version of THERMUS framework [2^ is used. A numerical 
procedure is applied to meet the above constraints simultaneously. Other choices of the freeze- 
out parameters will be discussed in the next section. The T, /i^, and 75 parameters used for 
different cm. energies are given in Table I. Here, some further details should be mentioned. We 
use quantum statistics, but disregard the non-zero widths of resonances. The thermodynamic 
limit for the calculations of the scaled variance is assumed, thus u reaches its limiting value. 
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and volume V is not a parameter of our model calculations. We also do not consider explicitly 
momentum conservation as it can be shown that it completely drops out in the thermodynamic 



limit 



171]. Excluded volume corrections due to a hadron hard core volume are not taken into 



account. They will be considered elsewhere 23|]. The standard THERMUS particle table 
includes all known particles and resonances up to a mass of about 2.5 GeV and their respective 
decay channels. Heavy resonances have not always well established decay channels. We re- 
scaled the branching ratios given in THERMUS to unity, where it was necessary to ensure a 
global charge conservation. Usually the resonance decays are considered in a successive manner, 
hence, each resonance decays into lighter ones until only stable particles are left. However, we 
need to implement another procedure when different branches are defined in a way that final 
states with only stable hadrons are counted. This distinction does not affect mean quantities, 
but for fluctuations it is crucial. To make a correspondence with NA49 data, both strong and 
electromagnetic decays should be taken into account, while weak decays should be omitted. 
Once a suitable set of thermodynamical parameters is determined for each collision energy. 
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[ MeV ] 


[ MeV ] 




GCE 
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MCE 
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2.32 


64.9 


800.8 


0.641 


1.025 


0.777 


0.578 


1 


020 


0.116 


0.086 


1.048 


0.403 


0.300 


4.86 


118.5 


562.2 


0.694 


1.058 


0.619 


0.368 


1 


196 


0.324 


0.192 


1.361 


0.850 


0.505 


6.27 


130.7 


482.4 


0.716 


1.069 


0.640 


0.346 


1 


203 


0.390 


0.211 


1.431 


0.969 


0.524 


7.62 


138.3 


424.6 


0.735 


1.078 


0.664 


0.334 


1 


200 


0.442 


0.222 


1.476 


1.060 


0.534 


8.77 


142.9 


385.4 


0.749 


1.084 


0.683 


0.328 


1 


197 


0.479 


0.230 


1.504 


1.126 


0.541 


12.3 


151.5 


300.1 


0.787 


1.097 


0.729 


0.320 


1 


185 


0.563 


0.247 


1.557 


1.271 


0.558 


17.3 


157.0 


228.6 


0.830 


1.108 


0.768 


0.318 


1 


174 


0.637 


0.263 


1.593 


1.393 


0.576 


62.4 


163.1 


72.7 


0.975 


1.127 


0.827 


0.316 


1 


147 


0.782 


0.298 


1.636 


1.609 


0.613 


130 


163.6 


36.1 


0.998 


1.131 


0.827 


0.313 


1 


141 


0.805 


0.305 


1.639 


1.631 


0.618 


200 


163.7 


23.4 


1.000 


1.133 


0.826 


0.312 


1 


140 


0.811 


0.307 


1.639 


1.636 


0.619 


5500 


163.8 


0.9 


1.000 


1.136 


0.820 


0.310 


1 


137 


0.820 


0.309 


1.640 


1.640 


0.619 



TABLE I: The chemical freeze-out parameters T, fis, 7s j and final state scaled variances in the GCE, 
CE, and MCE for central Pb-|-Pb (Au-|-Au) collisions at different cm. energies, ^/sJiN- 
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the scaled variance of negatively, positively, and all charged particles can be calculated using 
Eqs. (EE)- The Eqs. (IHlfTTl) lead to the scaled variance in the GCE, whereas Eqs. ( TT2lfT5l) 
correspond to the MCE (or CE) results. The cj^, u;^, uj'^^ in different ensembles are presented 
in Table I for different collision energies. The values of ^Jsnn quoted in Table I correspond 
to the beam energies at SIS {2A GeV), AGS (11. GeV), SPS (20A, 30A, 40A, 80A, and 
158A GeV), colliding energies at RHIC (^s^ = 62.4 GeV, 130 GeV, and 200 GeV), and LHC 
(^5^ = 5500 GeV). 

The mean multiplicities, (A^j), used for calculation of the scaled variance (see Eq. ([6])) are 
given by Eqs. (j2]) and ([3]) and remain the same in all three ensembles. The variances in Eq. ([6]) 
are calculated using the corresponding correlators (ANiANj) in the GCE, CE, and MCE. For 
the calculations of final state correlators the summation in Eq. (fT5l) should include all resonances 
R and R' which have particles of the species i and/or j in their decay channels. The resulting 
scaled variances are presented in Table I and shown in Figs. [TH as the functions of ^/sNN. 
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FIG. 1: The scaled variances for negatively charged particles, lo^ , both primordial and final, along the 
chemical freeze-out line for central Pb+Pb (Au+Au) collisions. Different lines present the GCE, CE, 
and MCE results. Symbols at the lines for final particles correspond to the specific collision energies 
pointed out in Table I. The arrows show the effect of resonance decays. 
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FIG. 2: The same as in Fig. [H but for 
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FIG. 3: The same as in Figs. [T]and[2l but for a; 
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At the chemical freeze-out of heavy-ion coUisions, the Bose effect for pions and resonance 
decays are important and thus (see also Ref. uj^^^ = 1.1, ujt^^ = 1.2, and u'^^ — 



^g.c.e. 



-5 "^g.c.e. 



^g.c.e. 
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1.4 1.6, at the SPS energies. Note that in the Boltzmann approximation and neglecting the 
resonance decay effect one finds cu^^.e. = ^t-c.e. = ^f^ce. = 1- 

Some quahtative features of the results should be mentioned. The effect of Bose and Fermi 
statistics is seen in primordial values in the GCE. At low temperatures most of charged hadrons 
are protons, and Fermi statistics dominates, i^^ce ' ^'tce < 1- other hand, in the limit 

of high temperature (low ^b/T) most charged hadrons are pions and the effect of Bose statistics 
dominates, cu^c.e.^ "^gx.e. > 1- Along the chemical freeze-out line, cu^c.e. is always slightly larger 
than 1, as TT^ mesons dominate at both low and high temperatures. The bump in uj^^ ^ for final 
state particles seen in Fig. [2] at the small collision energies is due to a correlated production 
of proton and vr"*" meson from A"'"^ decays. This single resonance contribution dominates in 
uj'^ce small collision energies (small temperatures), but becomes relatively unimportant at 
the high collision energies. 

A minimum in uj~^ for final particles is seen in Fig. [H This is due to two effects. As the 
number of negatively charged particles is relatively small, (A^-) <^ the low collision 

energies, both the CE suppression and the resonance decay effect are small. With increasing 
■s/snn, the CE effect alone leads to a decrease of uj~^, but the resonance decay effect only 
leads to an increase of uj~^. A combination of these two effects, the CE suppression and the 
resonance enhancement, leads to a minimum of uj~^. As expected, ujm.c.e. < ^c.e., as an energy 
conservation further suppresses the particle number fluctuations. A new feature of the MCE 
is the additional suppression of the fluctuations after resonance decays. This is discussed in 
Appendix A. 

IV. COMPARISON WITH NA49 DATA 
A. Centrality Selection 

The fluctuations in nucleus-nucleus collisions are studied on an event-by-event basis: a given 
quantity is measured for each collision and a distribution of this quantity is measured for a se- 
lected sample of these collisions. It has been found that the fluctuations in the number of 
nucleon participants give the dominant contribution to hadron multiplicity fluctuations. In the 
language of statistical models, fluctuations of the number of nucleon participants correspond 
to volume fluctuations caused by the variations in the collision geometry. Mean hadron multi- 
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plicities are proportional (in the large volume limit) to the volume, hence, volume fluctuations 
translate directly to the multiplicity fluctuations. Thus a comparison between data and predic- 
tions of statistical models should be performed for results which correspond to collisions with 
a flxed number of nucleon participants. 

Due to experimental limitations it is only possible to approximately measure the number 
of participants of the projectile nucleus, N^"-' , in flxed target experiments (e.g. NA49 at the 
CERN SPS). This is done in NA49 by measuring the energy deposited in a downstream Veto 
calorimeter. A large fraction of this energy is due to projectile spectators N^"-' . Using baryon 
number conservation for the projectile nucleus {A = N^"-' + N^"-') the number of projectile 
participants can be estimated. However, also a fraction of non-spectator particles, mostly 
protons and neutrons, contribute to the Veto energy pL4|. Furthermore, the total number 
of nucleons participating in the collision can fluctuate considerably even for collisions with 
a flxed number of projectile participants (see Ref. I2J]). This is due to fluctuations of the 
number of target participants. The consequences of the asymmetry in an event selection depend 
on the dynamics of nucleus- nucleus collisions (see Ref. 25| for details). Still, for the most 
central Pb+Pb collisions selected by the number of projectile participants an increase of the 
scaled variance can be estimated to be smaller than a few % 25| due to the target participant 
fluctuations. In the following our predictions will be compared with the preliminary NA49 data 
on the 1% most central Pb+Pb collisions at 20^-158^4 GeV [l^. The number of projectile 
participants for these collisions is estimated to be larger than 193. 



B. Modelling of Acceptance 

In the experimental study of nuclear collisions at high energies only a fraction of all produced 
particles is registered. Thus, the multiplicity distribution of the measured particles is expected 
to be different from the distribution of all produced particles. Let us consider the production 
of particles with the probability PiniN) in the full momentum space. If particle detection is 
uncorrelated, this means that the detection of one particle has no influence on the probability 
to detect another one, the binomial distribution can be used. For a flxed number of produced 
particles the multiplicity distribution of accepted particles reads: 

Paccin, N) = g"(l - g)^-" ■ f , (17) 

n! A — n)\ 
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where n < N and q is the probabihty of a single particle to be accepted (i.e. it is the ratio 
between mean multiplicity of accepted and all hadrons). Consequently one gets, n = q N , n^ — 



n 



q{l — q)N , where = "^^=0 Paccif^-, N) , for k = 1,2, . . . . The probability distribution 



P{n) of the accepted particles reads: 

oo 

P{n) = 5^P4.(iV) Pacc{n,N) . 

N=n 

The first two moments of the distribution P{n) are calculated as: 

oo N 



) = J2p,AN) J2^Paccin,N) = q-{N), (19) 

N=0 n=0 

oo N 

) = 5^P4.(iV) 5^n2p„,,(n,iV) = q'-{N') + q{l - q) ■ (N) , (20) 



in') = 



N=0 n=0 

where (fc = 1, 2, . . .) 



(iV') = X]iV^P4.(iV) . (21) 

N=0 

Finally, the scaled variance for the accepted particles can be obtained: 

CO ^ ^""'^ ~ ^""^^ = I - q +q-u;,^, (22) 
(n) 

where uj^t^ is the scaled variance of the P^-^iN) distribution. The limiting behavior of u agrees 
with the expectations. In the large acceptance limit (g ~ 1) the distribution of measured 
particles approaches the distribution in the full acceptance. For a very small acceptance (g ~ 
0) the measured distribution approaches the Poisson one independent of the shape of the 
distribution in the full acceptance. 

Model results on multiplicity fluctuations presented in Sec. Ill correspond to an ideal situ- 
ation when all final hadrons are accepted by a detector. For a comparison with experimental 
data a limited detector acceptance should be taken into account. Even if primordial particles at 
chemical freeze-out are only weakly correlated in momentum space this would no longer be valid 
for final state particles as resonance decays lead to momentum correlations for final hadrons. 
In general, in statistical models, the correlations in momentum space are caused by resonance 
decays, quantum statistics and the energy-momentum conservation law, which is implied in 
the MCE. In this paper we neglect these correlations and use Eqs. (fT8|) and (!22|) . This may 
be approximately valid for uj~^ and uj~ , as most decay channels only contain one positively (or 
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negatively) charged particle, but is certainly much worse for uj'^^, for instance due to decays of 
neutral resonances into two charged particles. In order to limit correlations caused by resonance 
decays, we focus on the results for negatively and positively charged hadrons. A discussion of 
the effect of resonance decays to the acceptance procedure and a comparison with the data for 
uj'^^ are presented in Appendix B. An improved modelling of the effect regarding the limited 
experimental acceptance will be a subject of a future study. 

C. Comparison with the NA49 Data for a; and uj^ 



Fig. m presents the scaled variances u and calculated with Eq. ( 122|) . The hadron- 
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FIG. 4: The scaled variances for negative (left) and positive (right) hadrons along the chemical freeze- 
out line for central Pb+Pb collisions at the SPS energies. The points show the preliminary data of 
NA49 



14l ]. Total (statistical+systematic) errors are indicated. The statistical model parameters T, 
Hb, and 75 at different SPS collision energies are presented in Table I. Lines show the GCE, CE, and 
MCE results calculated with the NA49 experimental acceptance according to Eq. (j22p . 



resonance gas calculations in the GCE, CE, and MCE shown in Figs. 1 and 2 are used for 
the cu^. The NA49 acceptance used for the fluctuation measurements is located in the for- 
ward hemisphere (1 < y^n) < ybeam, where y^ir) is the hadron rapidity calculated assum- 
ing pion mass and shifted to the collision cm. system [l^). The acceptance probabilities 
for positively and negatively charged hadrons are approximately equal, ~ q~ , and the 



numerical values at different SPS energies are: g"* 



0.038, 0.063, 0.085, 0.131, 0.163 at 
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^/sNN = 6.27, 7.62, 8.77, 12.3, 17.3 GeV, respectively. Eq. (12^ has the following property: 
if uj4t^ is smaller or larger than 1, the same inequality remains to be valid for uj at any value 
of < g < 1. Thus one has a strong qualitative difference between the predictions of the 
statistical model valid for any freeze-out conditions and experimental acceptances. The CE 
and MCE correspond to cu^.c.e. < ^t.e. < the GCE to ou^^.e. > 1- 

From Fig. H] it follows that the NA49 data for extracted from 1% of the most central 
Pb+Pb collisions at all SPS energies are best described by the results of the hadron-resonance 
gas model calculated within the MCE. The data reveal even stronger suppression of the particle 
number fluctuations. 



D. Dependence on the Freeze-out Parameters 



The relation E/N = 1 GeV 22| was used in our calculations to define the freeze-out con- 
ditions. It does not give the best fit of the multiplicity data at each specific energy. In this 
subsection we check the dependence of the statistical model results for the scaled variances on 
the choice of the freeze-out parameters. For this purpose we compare the results obtained for 
the parameters used in this paper (model A) with two other sets of the freeze-out parameters 
at SPS energies: model B 6|] and model C [7]. The corresponding values of T and are 
presented in Fig. 5. 
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0.211 
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0.334 
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0.232 


0.236 


12.3 


0.320 


0.318 


0.325 


0.247 


0.249 


0.248 


17.3 


0.318 


0.317 


0.321 


0.263 


0.264 


0.259 



TABLE II: Final state scaled variances calculated in the MCE for 4 vr acceptance using freeze-out 
conditions A, B, and C. 



The scaled variances c^m.c.e. c^m.c.e. calculated in the full phase space within the MCE 
vary by less than 1% when changing the parameter set. In the NA49 acceptance the difference 
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FIG. 5: Chemical freeze-out points in the T-fiB plane for central Pb+Pb collisions. The solid lines 
shows {E) I (N) = 1 GeV, the squares are from our parametrization (model A) and denote SPS beam 
energies from 20 A GeV (right) to 158^1 GeV (left), the full and open circles are the best fit parameters 



is almost completely washed out. The differences are somewhat stronger in the GCE and CE, 
but will not be considered here. 

E. Comparison of Distributions 

As discussed in Section II the multiplicity distribution in statistical models in the full phase 
space and in the large volume limit approaches a normal distribution. If the particle detection 
is modelled by the simple procedure presented in Section IV B then the results (11911221) are 
valid for any form of the full acceptance distribution P4^(A^). In the following we discuss the 
properties of the multiphcity distribution in the limited acceptance, P{n), (fTS!) and compare the 
statistical model results in different ensembles with data on negatively and positively charged 
hadrons. 

For the Poisson distribution in the full acceptance the summation in Eq.( |T8l) leads also to 
the Poisson distribution in the acceptance with the expectation value (n) = q{N): 



from reference 



(model B) and [7| (model C), respectively. 




n 



(23) 
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However, the same does not hold true for summation in Eq. (1181) being apphed to other forms 
of the distribution P{N). In particular, the normal distribution (fT6l) is transformed into the 
following: 

oo 

P{n) = Y,Pg{N) Pacc{n,N) , (24) 

N=n 

which is not anymore the Gauss one. It is enough to mention that a Gaussian is symmetric 
around its mean value, while the distribution (IHl) is not. 



The average number particles accepted by a detector is: 

oo 

(n) = ^nP{n) = q {N) = q p V , (25) 

n=0 

where p = {N)/V is the corresponding particle density. Hence, one can determine the volume 
to be 

q p 

In the following for each beam energy we adjust the volume to match the condition of Eq. (126!) 
for negatively (V^) and positively (V^) charged yields, separately. Note that values for the 
volume are about 10 — 20% larger than the ones in jfi, Q], which were obtained using a much 
less stringent centrality selection (here only the 1% most central data is analyzed). We find 
that the V~ and volume parameters deviate from each other by less than 10%. Deviations 
of a similar magnitude are observed between the data on hadron yield systematics and the 
hadron-resonance gas model fits. Here we are only interested in the shape of multiplicity distri- 
butions and do not attempt to optimize the model to fit simultaneously yields of positively and 
negatively charged particles. As typical examples the multiplicity distributions for negatively 
and positively charged hadrons produced in central Pb+Pb collisions at 40^4 GeV are shown 
in Fig. [6] at the SPS energy range. 

The bell-like shape of the measured spectra is well reproduced by the shape predicted by 
the statistical model. In the semi-logarithmic plot differences between the data and model 
lines obtained within different statistical ensembles are hardly visible. In order to allow for a 
detailed comparison of the distributions the ratio of the data and the model distributions to 
the Poisson one is presented in Fig. [71 

The results for negatively and positively charged hadrons at 20^4 GeV, 30A GeV, 40^4 GeV, 
80A GeV, and 158A GeV are shown separately. The convex shape of the data reflects the fact 
that the measured distribution is significantly narrower than the Poisson one. This suppression 
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FIG. 6: The multiplicity distributions for negatively (left) and positively (right) charged hadrons 



14| . The preliminary 



produced in central (1%) Pb+Pb collisions at AOA GeV in the NA49 acceptance 
experimental data (solid points) of NA49 [l^ are compared with the prediction of the hadron-resonance 
gas model obtained within different statistical ensembles, the GCE (dotted lines), the CE (dashed- 
dotted lines) and the MCE (sohd hues). 

of fluctuations is observed for botfi cfiarges, at all five SPS energies and it is consistent with 
the results for the scaled variance shown and discussed previously. The GCE hadron-resonance 
gas results are broader than the corresponding Poisson distribution. The ratio has a concave 
shape. An introduction of the quantum number conservation laws (the CE results) leads to the 
convex shape and significantly improves agreement with the data. Further improvement of the 
agreement is obtained by the additional introduction of the energy conservation law (the MCE 
results). The measured spectra surprisingly well agree with the MCE predictions. 



F. Discussion 



High resolution of the NA49 experimental data allows to distinguish between multiplicity 
fluctuations expected in hadron-resonance gas model for different statistical ensembles. The 
measured spectra clearly favor predictions of the micro-canonical ensemble. Much worse de- 
scription is obtained for the canonical ensemble and a strong disagreement is seen considering 
the grand canonical one. All calculations are performed in the thermodynamical limit which is 
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FIG. 7: The ratio of the multiphcity distributions to Poisson ones for negatively (upper panel) and 
positively (lower panel) charged hadrons produced in central (1%) Pb+Pb collisions at 20A GeV, 
30^ GeV, AO A GeV, 80A GeV, and 158^ GeV (from left to right) in the NA49 acceptance Q. The 
preliminary experimental data (solid points) of NA49 [IJ] are compared with the prediction of the 
hadron-resonance gas model obtained within different statistical ensembles, the GCE (dotted lines), 
the CE (dashed-dotted lines), and the MCE (solid lines). 

a proper approximation for the considered reactions. Thus these results should be treated as 

n 

a first observation of the recently predicted [15| suppression of multiplicity fiuctuations due to 
conservation laws in relativistic gases in the large volume limit. 

A validity of the micro-canonical description is surprising even within the framework of the 
statistical hadronization model used in this paper. This is because in the calculations the pa- 
rameters of the model (e.g. energy, volume, temperature and chemical potential) were assumed 
to be the same in all collisions. On the other hand, significant event-by-event fiuctuations of 
these parameters may be expected. For instance, only a part of the total energy is available for 
the hadronization process. This part should be used in the hadron-resonance gas calculations 
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while the remaining energy is contained in the collective motion of matter. The ratio between 
the hadronization and collective energies may vary from collision to collision and consequently 
increase the multiplicity fluctuations. 

The agreement between the data and the MCE predictions is even more surprising when 
the processes which are beyond the statistical hadronization model are considered. Examples 
of these are jet and mini-jet production, heavy cluster formation, effects related to the phase 
transition or instabilities of the quark-gluon plasma. Naively all of them are expected to 
increase multiplicity fluctuations and thus lead to a disagreement between the data and the 
MCE predictions. A comparison of the data with the models which include these processes is 
obviously needed for significant conclusions. Here we consider only one example. 



In Ref. HI] a non-monotonic dependence of the relative fluctuations. 



has been suggested as a signal for the onset of deconfinement. Here S and E denote the system 
entropy and thermalized energy at the early stage of collisions, respectively. This prediction 
assumes event-by-event fluctuations of the thermalized energy, which results in the fluctuations 
of the produced entropy. The ratio of the entropy to energy fluctuations ( 1271) depends on the 
equation of state and thus on the form of created matter. The is approximately independent 
of collision energy and equals about 0.6 in pure hadron or quark-gluon plasma phases. An 



increase of the Re ratio up to its maximum value. Re ~ 0.8, is expected 



ll| in the transition 



domain. Anomalies in energy dependence of the hadron production properties measured in 



central Pb+Pb collisions [26] indicate 271] that the transition domain is located at the low 



CERN SPS energies, from 30^4 to SOA GeV. Thus an anomaly in the energy dependence of 



multiplicity fluctuations is expected in the same energy domain 

In any case the fluctuations of the thermalized energy will lead to additional multiplicity 
fluctuations ("dynamical fluctuations"). The resulting contribution to the scaled variance can 
be calculated to be: 

^-,yn = Re (^-) ^ • (28) 

The above assumes that the mean particle multiplicity is proportional to the early stage en- 
tropy. In order to perform a quantitative estimate of the effect the fluctuations of the energy 

n n 

of produced particles were calculated within the HSD ^] and UrQMD [29[ string-hadronic 
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models. For central (impact parameter zero) Pb+Pb collisions in the energy range from 30A 



to 80A GeV we have obtained, a/ {5EY / E < 0.03. The number of accepted negatively charged 
particles is (jiJ) ^ 30 at 40^4 GeV (see Fig. [7]). Thus, an increase of the uj due to the "dynam- 
ical fluctuations" estimated by Eq. (1281) is oJ'^y^ < 0.02 for R^. = 0.6, and it is smaller than the 
experimental error of the preliminary NA49 data of about 0.05 |l4|. In particular, an additional 
increase due to the phase transition, Au^y^ ~ 0.005, for Re = 0.8, can be hardly observed. 

In conclusion, the predicted HI] increase of the scaled variance of the multiplicity distribution 
due to the onset of deconfinement is too small to be observed in the current data. These data 
neither confirm nor refute the interpretation [27] of the measured {2^ anomalies in the energy 
depedence of other hadron production properties as due to the onset of deconfinement at the 
CERN SPS energies. 

More differential data on multiplicity fluctuations and correlations are required for further 
tests of the validity of the statistical models and observation of possible signals of the phase 
transitions. The experimental resolution in a measurement of the enhanced fluctuations due 
to the onset of deconfinement can be increased by increasing acceptance. For example, uj^y^ oc 
oc q. The present aceptance of NA49 at AOA GeV is about q = 0.06 and it can be increased 
up to about q = 0.5 in the future studies. This will give a chance to observe, for example, the 
dynamical fluctuations discussed in Ref. jllj]. The observation of the MCE suppression effects 
of the multiplicity fluctuations by NA49 was possible only because a selection of a sample 
of collisions without projectile spectators. This selection seems to be possible only in the 
fixed target experiments. In the collider kinematics nuclear fragments which follow the beam 
direction can not be measured. 

On the model side a further study is needed to improve description of the effect of the 
limited experimental acceptance. Further on, a finite volume of hadrons is expected to lead to 
a reduction of the particle number fluctuations [23.]. A quantitative estimate of this effect is 
needed. 



V. SUMMARY 

The hadron multiplicity fluctuations in relativistic nucleus-nucleus collisions have been pre- 
dicted in the statistical hadron-resonance gas model within the grand canonical, canonical, and 
micro-canonical ensembles in the thermodjTiamical limit. The microscopic correlator method 
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has been extended to include three conserved charges - baryon number, electric charge, and 
strangeness - in the canonical ensemble, and additionally an energy conservation in the micro- 
canonical ensemble. The analytical formulas are used for the resonance decay contributions 
to the correlations and fluctuations. The scaled variances of negatively, positively, and all 
charged particles for primordial and final state hadrons have been calculated at the chemical 
freeze-out in central Pb-|-Pb (Au-|-Au) collisions for different coUision energies from SIS to LHC. 
A comparison of the multiplicity distributions and the scaled variances with the preliminary 
NA49 data on Pb-|-Pb coUisions at the SPS energies has been done for the samples of about 
1% of most central collisions selected by the number of projectile participants. This selection 
allows to eliminate effect of fluctuations of the number of nucleon participants. The effect of 
the limited experimental acceptance was taken into account by use of the uncorrelated particle 
approximation. 

The measured multiplicity distributions are significantly narrower than the Poisson one and 
allow to distinguish between model results derived within different statistical ensembles. The 
data surprisingly well agree with the expectations for the micro-canonical ensemble and exclude 
the canonical and grand-canonical ensembles. Thus this is a first experimental observation of 
the predicted suppression of the multiplicity fluctuations in relativistic gases in the thermody- 
namical limit due to conservation laws. 
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APPENDIX A: RESONANCE DECAYS IN THE MCE 



A comparison of the primordial scaled variances with those for final hadrons demonstrates 
that the fluctuations generally increase after resonance decays in the GCE and CE (see more 
details in Ref. |l6!]), but they decrease in the MCE. In order to understand this effect let us 
consider a toy model (tt^, tt^, p'^)-system with a zero net charge, Q = 0. Due to this last 
condition there is a full symmetry between positively and negatively charged pions, and thus 
LV~^ = uj~ . Each p°-meson decays into a Tr+vr^-pair with 100% probability, i.e. = 1 and 
(n_)pO = 1. The predictions of the CE and MCE for (vr^, vr", p°)-system are shown in Fig. [HI 
One observes that p° decays lead to an enhancement of in the CE, and to its suppression 




60 120 180 60 120 180 

T (MeV) T (MeV) 



FIG. 8: The scaled variance for negatively charged particles, lo^ , in the toy (vr^, vr^, /9'^)-system 
with Q = in the CE (left) and MCE (right) as a function of temperature. The temperature 
interval corresponds to that presented in Table I. The dotted lines correspond to primordial vr^-meson 
fluctuations, the solid lines - to those after decays. At small T the contribution to the pion 
number fluctuations is negligible, and it increases with T. The contribution from p^ decays to lu^ is 
positive in the CE and negative in the MCE (see the text for details). 

in the MCE. In the CE one finds from Eqs. ([2]) and (ITSl) for (tt^, tt^, p°)-system: 

(iv_) = (iv;_) + (iv,o), ((Aiv_)2)_ = ((Aiv;.)^)_ + ((Aiv,o)^)_ . (ai) 
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Note that the average multiphcities, {N*_) and (Npo), remain the same in the CE, and the 
MCE. From Eq. (jSH) it follows: 



(iV_) 



(A2) 



The uj~* is essentially smaller than 1 due to the strong CE suppression (see Fig. [HI left). On 
the other hand, there is no CE suppression for p° fluctuations, uf^ = ujfce — 1- Therefore, 
one finds that ^^c°e./^c.e. > 1; from Eq. (]A2p it immediately follows, uj~^^ > cj^* . Note that 
a;~*g = w^^c g = 1, thus there is no enhancement of uj~^ g due to decays. In the MCE the 
multiplicity (A^-) remains the same as in Eq. ( lAll) . The variance {{AN^)'^)m.c.e. is, however, 
modified because of the anti-correlation between primordial 7r~* and p° mesons in the MCE. 
From Eq. (fTSll one finds for our (tt^, tt^, p°)-system. 



{{AN_y)m.c.e. = {{AN:.y)m.c.e. + {{AN pof) + 2 (AA^^ ANpo)^,,,,, . (A3) 

The last term in Eq. ( ]A3I) appears due to energy conservation in the MCE (it is absent in the 
CE). This term is evidently negative, which means that an anti-correlation occurs. A large 
(small) number of primordial pions, AN*^ > (< 0), requires a small (large) number of p"- 
mesons, AN^o < (> 0), to keep the total energy fixed. Anti-correlation between primordial 
pions and p°-mesons makes the n~ number fluctuations smaller after resonance decays, i.e. 
'^m.c.e. < ^rn*c.e.^ depicted in Fig. M (right). The same mechanism works in the MCE for the 
full hadron-resonance gas. 



APPENDIX B: ACCEPTANCE EFFECT FOR ALL CHARGED PARTICLES 

In order to better understand an influence of the momentum correlation due to resonance 
decays on the multiplicity fluctuations we define a toy model. Let us assume that there are 
two kinds of particles produced. The first kind (A^) is either stable or originates from decay 
channels which contain only one particle of the type we are set to investigate, while the second 
kind (M) produces 2 particles of the selected type. This is described by the (unknown) proba- 
bility distribution P^t^ {N, M). We further assume that for both types of particles, A^ and M, 
separately the acceptance procedure defined by Eq. (ITTl) is applicable. We also assume that 
once particle M is inside the experimental acceptance, both its decay products will be so as 
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well. Hence, the average number of observed particles will be: 

oo oo N M 



N=0 M=0 

This leads immediately to: 



n=0 m=0 



{n) = q- {N) + 2 (M) 



(Bl) 



fB2) 



One finds the second moment, 

oo oo 

') = P^AN,M) (n + ^mf Pace in,N) Pace {m,M) . 



N M 



in 



(B3) 



N=0 M=0 n=0 m=0 

Making use of the relation fl20|) one obtains: 



(n^) = q(l-q).{N) + ■ {N^) + V • {NM) + 4g (1 - g) ■ (M) + V " (M^) • (B4) 



Thus, for the scaled variance it follows: 



{n) 



q ■ 



{N) + 4 (M) 
{N) + 2 (M) 



(B5) 



where is obtained from the case q = 1 and corresponds to the original distribution 
P4t^{N,M). For the second limiting case of Eq. (IBSp . g — 0, one finds a scaled variance 
which corresponds to that of two uncorrelated Poisson distributions with means q{N) and 
q{M), respectively. In this case, all primordial correlations due to energy and charge conser- 
vation or Rose (Fermi) statistics are lost, but particles produced by resonances of type M are 
still detected in pairs. In general case, by k we denote the fraction of particles originating from 
decays (always 2 relevant daughters) of particle kind M, hence. 



{N) = {1-k) {Ntot) , 
Finally, one finds for the scaled variance: 



(M) 



^ (Ntot) . 



(B6) 



q-LU^^ + (l-g)(l + fc) . 



(B7) 



From the hadron-resonance gas model we can estimate the fraction k of the final yield which 
originates from decays of resonances into 2 (or more) charged particles. From Fig. [9] (left) we 
find the fraction of the charged particle yield k to be from 35% (20 AGeV) to 45% (158 AGeV) 
in the SPS energy range (and about 10% for positively and negatively charged particles). 
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FIG. 9: Left. Fraction of total yield originating from resonance decays with 2 or more positively('+'), 
negatively('-'), or all charged('ch'), particles. Parameters are taken from parametrization in section 
III. Right. The scaled variances for all charged hadrons along the chemical freeze-out line for central 
Pb+Pb collisions at the SPS energies. The points show the preliminary data of NA49. The statistical 
model parameters T, fis, and 75 at different SPS collision energies are presented in Table I. The lines 
show the MCE results calculated with the NA49 experimental acceptance according to Eq. (j22p (lower 
line) and Eq. ()B7p ( upper line). 

For the definition of decay channels see section III. Examples of two-particle decay channels: 
TT^ + TT^, would be counted as two particle decay in 'ch', but neither in '+' nor in 
—>■ p+7i^, would contribute to 'ch' and '+', but not to The assumption that both decay 
products are detected is certainly not justifiable for small values of the total acceptance q, hence 
Eq. (IB7I) overestimates the effect. However, this consideration will give a useful upper bound 
(see Fig. [9], right). The typical width of decays is comparable to the width of the acceptance 
window, therefore, about half of all decays will leave one (or both) decay product missing. Yet 
the same 50% will be contributed from decays whose parents are outside the acceptance but 
contribute to the final yield. Hence, one expects no change in average multiplicity, but a sizable 
effect on fluctuations. 
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